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We present the calculation of the coherent spectral functions and density of states (DOS) for 
excitonic systems in the frame of the three dimensional extended Falicov-Kimball model. By using 
gauge-invariant U(l) transformation to the usual fermions, we represent the electron operator as 
a fermion attached to the U(l) phase-flux tube. The emergent bosonic gauge field, related to the 
phase variables is crucial for the Bose-Einstein condensation (BEC) of excitons. Employing the 
path-integral formalism, we manipulate the bosonic and fermionic degrees of freedom to obtain the 
effective actions related to fermionic and bosonic sectors. Considering the normal and anomalous 
excitonic Green functions, we calculate the spectral functions, which have the forms of convolutions 
in the reciprocal space between bosonic and fermionic counterparts. For the fermionic incoherent 
part of the DOS we have found the strong evidence of the hybridization-gap in DOS spectra. 
Furthermore, considering Bogoliubov coherence mechanism, we calculate the coherent DOS spectra. 

For the coherent normal fermionic DOS, there is no hybridization-gap found in the system due to 
strong coherence effets and phase stiffness. The similar behavior is observed also for the condensate 
part of the anomalous excitonic DOS spectra. We show that for small values of the Coulomb 
interaction, fermionic DOS exhibits a Bardeen-Cooper-Schriffer (BCS) -like double-peak structure. 

In the BEC region of the BCS-BEC crossover, the double-peak structure disappears totally for both: 
coherent and incoherent DOS spectra. We discuss also, temperature dependence of DOS functions. 


I. INTRODUCTION 


The attractive Coulomb interaction between conduc¬ 
tion band electrons and valence band holes plays a cru¬ 
cial role in forming the bosonic electron-hole (e-h) bound 
pairsi (excitons). The Bose-Einstein condensation of the 
excitons represents a very interesting subject in solid 
state physics. Mostly, the e-h systems are realized in pho- 
toexcited semiconducting materials and the properties of 
such systems strongly depend on the e-h density, tem¬ 
perature of the system, pressure, and other important 
physical quantities. There is a number of works, both 
experimental and theoretical, where all those effects are 
investigated intensively— — Low density system of exci¬ 
tons behaves like the usual Bose gas, and, at cryogenic 
temperatures, the BEC transition can be expected— £ 
On the other hand, the high density system of bound 
e-h pairs behaves like the system of usual Cooper pairs 
in superconductors,^ In this limit we have the Bardeen- 
Cooper-Schrieffer (BCS) state of e-h pairs— £ Despite 
many experimental investigations to observe the coher¬ 
ent exciton condensates^— there is not yet definitive 
evidence for such states. Therefore, an expected BCS- 
BEC crossover represents actually a fascinating problem 
typical to the excitonic systems—Especially, it is in¬ 
teresting from the viewpoint of the difference from the 
similar crossover in superconductors or trapped atomic 
Fermi gases—The transition to the e-h pair condensed 
phase, in the limit of weak-coupling, is related to the 
relative motion between electrons and holes*^ implying 
the BCS regime and is in contrast to the case of strong¬ 
coupling regime, when the BEC state is related to the 
motion of the center of mass of excitons— In the whole 
BCS-BEC transition region, the e-h mass difference leads 


to a large suppression of the BEC transition tempera¬ 
ture, which is proved to not be same as the temperature 
of excitonic pair (EP) formation, and hence the excitonic 
insulator (El) phase— Recently, it is shown theoretically 
that the excitonic insulator and the excitonic condensate 
are not exactly the same states of the matter— The 
author in Refs. [24j and [2^ shows from general considera¬ 
tions that in the low density limit of the excitonic pairs, 
the critical temperature of excitonic condensation should 
be much smaller than the temperature of EP formation, 
in contrast to the previous treatments*^ - — where the El 
state is ad hoc associated with the BEC state of excitons 
as the same. Similarly, in Ref. [l^ it is shown that the El 
state is an excitonium state, where incoherent e-h bound 
pairs are formed and, furthermore, at the lower tem¬ 
peratures, the BEC of excitons appears in consequence 
of the reconfiguration and coherent condensation of the 
preformed excitonic pairs. Obviously, in the low den¬ 
sity limit, the gas of free excitons undergoes the BEC 
phase transition at the very low temperatures, and in 
general the BEC temperature transition line is not coin¬ 
ciding with that of EP formation. The Bose condensation 
of excitonic pairs is possible only when the macroscopic 
phase coherence is attended by the system. 2 - However, 
the experimental evidence of the existence of two distinct 
transition temperatures, for the general case of a three 
dimensional (3D) bulk system, is yet lacking in the liter¬ 
ature. In fact, the question, whether a true, “coherent” 
BEC transition is present in the system of excitons, (in 
the case of high exciton density limit) is still ambiguous. 
The experimental proof of it, is a very cumbersome prob¬ 
lem, because of dominant role of quantum fluctuations at 
low temperatures k B T < E^ y (with E e ^ y being the bind¬ 
ing energy of a Mott-Wannier exciton), when very large 
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zero-point oscillations are present. 

The continuing growing interest, to the problem of the 
coherent excitonic condensates motivates us to calculate 
excitonic spectral functions and density of states, as a 
direct probe mechanism, to compare the results with the 
high-resolution studies of angle-resolved photoemission 
spectroscopy (ARPES)2£~—, and with spectral weight 
measurements on the excitonic materials at the very low 
temperatures. 

The general strategy of our calculations is based on the 
effective actions method. First, we transform the initial 
total action of the system to a gauge invariant form, by 
applying the U(l) gauge transformation to the fermion 
operators. As a result, the electron appears in the the¬ 
ory like a composite object of that of the fermion with 
the attached U(l) phase-“flux-tube”. The electron fac¬ 
torization in terms of two variables has an unprecedented 
impact on the whole theory. Then we integrate out the 
phase variables and we get the effective fermionic action 
in the theory and we derive a set of self-consistent equa¬ 
tions for the El state. Furthermore, we discuss shortly 
the results of the quantum rotor model obtained after the 
integration of fermionic degrees of freedom. 

The path-integral formalism elaborated here permits 
to calculate the correlation functions in the system and, 
as a result, we obtain the expressions of normal (both in¬ 
coherent and coherent) and anomalous excitonic spectral 
functions, and shapes of density of states (DOS) corre¬ 
sponding. A special attention is payed, when calculating 
the bosonic phase-stiffness DOS function, which is nega¬ 
tive. Furthermore, it is crucial for calculation of the co¬ 
herent normal and anomalous excitonic DOS functions. 
Namely, considering the expressions of the bosonic and 
fermionic DOS functions, we calculate the total, phase 
coherent DOS functions, as convolutions from bosonic 
and fermionic counterparts. 

For the incoherent partial normal fermionic DOS func¬ 
tions at T = 0 we obtain a hybridization-gap in the exci¬ 
tation spectra as a direct consequence of the presence of 
the Hartree-type gap in the single-particle energy scales. 
For the case of coherent normal and anomalous excitonic 
DOS functions, this hybridization-gap is absent totally 
for all frequency modes and for all values of the Coulomb 
interaction. This is due of strong coherence effects in 
the strongly correlated fermion system at low tempera¬ 
tures and at low densities. For the anomalous excitonic 
DOS function, we found that the hybridization-gap is ab¬ 
sent for the case of small and intermediate values of the 
Coulomb interaction parameter, but there is a finite con¬ 
stant small gap that is opening in the strong coupling 
regime, signaling the passage to the SC (BEC) side of 
the SC-SM phase transition (the BCS-BEC crossover). 

The paper is organized as follows: The paper is 
organized as followsdn Section |TT] we introduce the 
model Hamiltonian and we discuss the main calculation 
schemes. In Section uni we describe the El state, with a 
short discussion about important energy scales in the sys¬ 
tem. In Section hyi we present the analytical calculation 


of spectral functions and DOS functions of the system 
in consideration. At the end of the Section na we show 
the results of numerical evaluations for calculated DOS 
functions. 


II. THE METHOD 
A. The EFKM Hamiltonian 

As the model for study the excitonic condensation 
at low temperatures we have chosen the two-band ex¬ 
tended Falicov-Kimball model (EFKM) ) 16 ' 30- — due to 
its large applicability for treatment of the electronic 
correlations ; 38 ' 39 The Hamiltonian of the EFKM model 
is 

H = —t c ^ [c(r)c(r') + h.c] — p ^ n(r) — 

<r,r') r 

-tf + h - c -\ + e ° 2 €J ™( r ) + 

(r,r'> r 

+U Y,\ [ n2 ( r ) “ ft2 ( r )] • (!) 

r 

Here /( r) (c(r)) is the / (c) electron creation operator at 
the lattice position r, the summation (r, r') is over near¬ 
est neighbors (n.n) sites on the 3D cubic lattice. The 
short hand notations are introduced n{ r) = n c (r) + n/(r) 
and h{ r) = n c (r) — n/(r) in order to simplify the calcu¬ 
lations. Next, t c is the hopping amplitude for c-electrons 
and e c is the corresponding on-site energy level parame¬ 
ter. Similarly, tf is the hopping amplitude for /-electrons 
and £f is the on-site energy level parameter for /-orbital. 
For t c tf <0 ( t c tf > 0) we have a direct (indirect) band 
gap semiconductor. The on-site (local) Coulomb interac¬ 
tion U in the last term of the Hamiltonian in Eq. JT]) plays 
the coupling role between two bands. As we will see later 
on, the strength of the local Coulomb interaction will 
tune the semi-metal (SM)- semiconductor (SC) transi¬ 
tion in the system. The chemical potential p is p = /i — e, 
where e = (e c + e/) /2. Naturally, we adjust the chemi¬ 
cal potentials fj ,/ and p, c in order to maintain the number 
of electrons in / and c orbitals separately. Then, the 
equilibrium value of the chemical potential n = /if = /jl c 
in Eq.([l]) will be determined from the half-filling con¬ 
dition, i.e., we suppose that ( n c (r )) + (n/(r)) = 1. We 
will use t c = 1 as the unit of energy and we fix the band 
parameter values e 0 = 0 and tf = —1. The Fermi en¬ 
ergy level is assumed to be situated at the level of the 
c -band, thus tp = 0. For the / -band hopping ampli¬ 
tude tf we consider the values tf = —0.3 and tf = —0.1 
corresponding to the heavy hole and light hole / -bands. 
Throughout the paper, we set &b = 1, = 1 and lat¬ 

tice constant a = 1. Also, we keep the frequency symbol 
v for fermions and w - for bosons, throughout the pa¬ 
per. In the case of degenerated / and c -bands, i.e., 
when tf = e c and tf = t c , the EFKM model reduces to 
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the standard Hubbard models The principal advantage 
of the EFKM, in comparison with the genuine Falicov- 
Kimball model (FKM), 3 - is that it is taking into account 
the direct nearest-neighbors / electron hopping a ^ ^ . (tf) 
and it can be shown2& that the El state is unstable when 
the pure FKM is approached. 

In fact, the EFKM Hamiltonian in Eq.© is equivalent 
to the asymmetric Hubbard model, if we associate to 
the orbitals c and / the spin variables, by replacing the 
fermionic Hilbert space with the pseudo-fermionic one 
and then linearizing the interaction term via the bosonic 
states (see in Ref. [ 13 ) ■ 


1. The partition function 


the exponent of the partition function in Eq. d2j , we have 


exp 




dr 



—iV (rr)n(rr) 


( 5 ) 


Furthermore, we combine the exponent in Eq.© with 
the effective chemical potential term, linear in total elec¬ 
tron density n(r) (see the second term, in Eq. ([]])). Then, 
we decompose the variables V{r r) into a static and the 
periodic part 


V (rr) = Vo (r) + V(rr), (6) 


The Hamiltonian in Eq. m is containing two sepa¬ 
rate quadratic terms and is suitable for decoupling by 
functional path integration methodj^tLS We employ the 
imaginary-time, fermionic path-integral method and we 
introduce the fermionic Grassmann variables^ / (rr) and 
c(rr) at each site r and each time r varying in the interval 
0 < t < j3, where /3 = 1/T with T being the thermo¬ 
dynamic temperature* 43 The time-dependent variables 
/(rr) and c(rr) are satisfying the anti-periodic boundary 
conditions x(r r) = — x[vt + (3), where x = f or c. 

The grand canonical partition function of system of 
fermions written as a functional integral over the Grass¬ 
mann fields is 


where, from time-periodicity of V(rr), it follows the re¬ 
lation Jq drF(rr) =0. As a result, the integration over- 
variables V(r t) becomes the integration over the scalar- 
static variables Vo(r) and the integration over the peri¬ 
odic field V’(rr): 


/ 


[DV] 



DV 


( 7 ) 


For the periodic part in Eq.©, we introduce U(l) phase 
field variables </>(rr) by using a Faraday-type relational 


V(vt) 


d<j>(rT) 

dr 


0(r t ). 


( 8 ) 


Zgc = J [DcDc] [DfDf] (2) 

where the action in the exponent is given in the path- 
integral formulation in the form 

x r P 

S[c,c,f,f]= ^2 Sb[x,x\+ / dr H{t). (3) 


*=f,c 

Here, Sn [/, /] and Sb[c, c] are fermionic Berry terms. 
They are defined as follows 


rP 

5_b[x, a:] = / dra;(rr)i:(rr) 

^ Jo 


( 4 ) 


Thus, for the dynamic part, we transform the integration 
over the gauge variables V(r r), into an integration over- 
generic phase variables /(rr) 


/ 


VV 


[Vcf] .... 


( 9 ) 


The periodicity of V (rr) implies that </>(r/3) = </>(rO). 
The integration measure, in Eq.©, over the variables 1 f> 
is defined as 



( 10 ) 


and i(rr) = d T x(r r) is the time-derivative . 


2. Decoupling of term proportional to n 2 

We will decouple the quadratic density terms in 
Eq.© using the Hubbard-Stratonovich linearisation 
procedure!- and by introducing new variables V (rr) and 
p(rr), conjugated respectively to the density terms n(rr) 
and h(rr). For quadratic term, proportional to n 2 (rr) in 


where the notations 4>i and </>/ mean the initial and fi¬ 
nal paths. The path-integral in Eq. m could be trans¬ 
formed into path integration over the compact U(l) 
group-manifold, since the electromagnetic group U(l), 
governing the phase field, is compact, i.e., />(rr) has the 
topology of a circle (S/), thus we have a non-homotopic 
mapping of the configuration space onto the U(l) gauge 
group Si —> 17(1). The paths, which loop around a circle 
in different number of times, are in different homotopy 
classes and they cannot be continuously deformed into 
one another. All these paths can be characterized by 
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their proper winding numbers^ m{ r). Any two paths 
which have different winding numbers, cannot be con¬ 
tinuously transformed one to another, and in order to 
include all possible phase path contributions, we have to 
sum over all topologically inequivalent phase configura¬ 
tions described by their winding numbers.— Accordingly, 
the path-integral in Eq. dlOl) is transformed as 

J m... = j pty]..., (li) 

where the integration measure is now 

[ [D(p]...= Y I 

l‘cp(r/3)=(p 0 (r)+2Trm(r) 

x / JJd^(rr)... . (12) 

J <p(rO)=ipo(r) r 

In performing the integration over the phase field one 
should take into account that the field configurations sat¬ 
isfy the boundary conditions^ 

tp(r/3) — y>(rO) = 2nm(r). (13) 

Thus, integration over all phases 4>{yt) amounts the in¬ 
tegration over the /3-periodic field </?(rr) and the sum¬ 
mation over a set of U(l) integer winding numbers m(r). 
For the scalar static part Fo(r), we get the following func¬ 
tional integral 

[DV 0 \ dr ~ Ys £ 1+iVt J ( r )[ n (r r )-v] . (14) 


The saddle-point value of Vo(r) is given by 



(15) 


3. Decoupling of term proportional to h 2 


The decoupling of the quadratic term proportional to 
n 2 (rr) in the exponent of the partition function in Eq. @ 
is also straightforward. We obtain 



After combining the expression in the exponent in 
Eq. (1181) with the similar linear term in the expression of 
the Hamiltonian in Eq. Jl|) (see the forth-term in Eq. (JT])) , 
we have 


[Dg\e 


J2 r fo dT ~ g2 4 ( U T) +e (rr) [n(rr)— ‘ c 2t f f ] 


(19) 


The saddle-point evaluation for g gives 



( 20 ) 


where n = (n(rr)) is the average of the particle den¬ 
sity difference function. As a result of the decoupling, 
we obtain a “Zeeman”- like contribution to the partition 
function 


exp 


~Y / dTiLnn( rr) 


( 21 ) 


with the attached effective chemical potential = 

e c~£f Un 

2 2 • 


4- Linearized, action with phase-field contribution 


where n is total average particle density n = n c + rif 
(furthermore, we will fix n as equal to 1, corresponding 
to the case of half-filling). 

Thereby, after decoupling of the quadratic term pro¬ 
portional to rr in the Hamiltonian in Eq. we get a 
contribution to the partition function in Eq.([2]) in the 
form 


exp 


5 [<p]-Y I d,Tg n n{YT) 


(16) 


To summarize, the grand canonical partition function 
of the system, after of both procedures of decoupling, is 

Zgg = J [ DcDc} [DfDf] [Dip] e~ s ^° (22) 

where the action <S[c, p\ in the exponent is given 

by 

«s[c, C, /, f,<p]=S [<p] + Y Sb x ] 

x=f,c 


where the emergent phase-only action <S[y>] is given as 



iu v<rT) 


i(p(rr)n(rT ) 


(17) 


and the effective chemical potential g n , attached to the 
total density operator in Eq. dlGI) . is given in the form 
Mn = — M- 



dr [c(rr)c(r'r) + h.c.] 



dr [/(rr)/(r'r) + h.c] 



dr [g n n(rT ) + ^n(rr)] . 


( 23 ) 
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After the Hubbard-Stratanovich linearisation, we got 
the total action of the system that is linear in terms 
of fermion densities and contains in addition a phase- 
dependent term £[</?], and also the terms, proportional 
to the effective chemical potentials \i n and pa. 


B. The U(l) transformation 

In the perspective to treat the local and non-local 
correlations in our excitonic system, it is important to 
separate the U(l) gauge degrees of freedom related to 
the phase sector. To this end, we perform the local 
gauge transformation to new fermion Grassmann vari¬ 
ables /(rr) and c(rr). Meanwhile, this procedure will au¬ 
tomatically eliminates also the last imaginary term, ap¬ 
pearing in the expression of the phase action in Eq. (El- 
For the electrons of / and c-orbitals, the U(l) trans¬ 
formation is 



FIG. 1: (Color online) Functional integration procedure. 


x(r t) 

= u(p)- 

x(yt) 

x(r t) 


x(r r) 


(24) 


where U(<p) is the U(l) transformation matrix lA(p) = 
I ■ cos p{vr) + ia z ■ sin p(yt) with the unit matrix I and 
a z being the Pauli matrix. The variables x = /, c, and we 
used the bosonic phase variables p introduced in Eqs. m 
and m ■ After transformations given in Eg. EH) , we 
obtain the total action of the system in the U(l) gauge 
invariant form (for comparison, see the action in Ea. (123l) 
before transformations) 


From this form of the partition function we will generate 
the effective actions for fermions and for bosonic phase 
sector (see the general procedure presented in Fig. [[]). 


III. EFFECTIVE ACTION FOR FERMIONS 


By following left-lowest root, presented in Fig. ??, we 
will integrate out the phase variables. We obtain 


<S[c, c, /, /, p\ = S 0 [< p } + ^2 Sb [x, x] 

X=/>C 


Z= [: DcDc] 


VfVf 


-Seff [c,C,/,/] 


(28) 


—t / dl ~ c(rr)c(r'r)e z M rr) r )] + h.c. 

/ / \ J 0 
<r,r') 


/ dr 


/(rr)/> , T)e- i M rT) - v(r ' T) ] + h.c. 
rP 

+ y2 dr [Unnirr) + ^n(rr)] . 

„ Jo 


(25) 


Now, t and t in Eci. (l25l) are respectively / -band and c 
-band fermion transfer integrals. We got in Ea. (l25l) also, 
a new, emergent, quadratic phase action iSqM 


5 oM = 


dr 


P 2 {vt) 2/i . ' 

—— lu‘ f{ ’ T) 


(26) 


The partition function of the system in the new variables 
/ and c is 


Z GG = / [D'cDc] \DfDf] [Dp] e-sfxJJM 


(27) 


where the effective, phase-averaged fermionic action in 
the exponent is 


«S eff [d, c, /, /] = -1 nj [Dp\e- S & a 'f’?^. (29) 


The Fourier transformation of fermionic variables / and 
c is 


x(rr) = ^j2 x ^y ikr ~ unT) ( 3 °) 

k,is n 


where N , is the number of lattice sites, and v n = 
7 t(2 n + l)/(3 are the Fermi-Matsubara frequencies^ with 
n = 0, Al, ±2,... . 

Then, the effective phase-averaged fermionic action of 
the system in the Fourier space takes the following form 


.Seff 


’ ^ = syE [ 5 kOn), /kW 

k Vn 


xQ x (k, v n ) 


CkOn) 

/k(t'n) 


(31) 
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Here, Q 4 (k, u n ) is the inverse of the Green function ma¬ 
trix 


S _1 (k,i/ n ) 


El{Vn) -A \ 

-A Ei(u n ))’ 


(32) 


where we have introduced the local excitonic order pa¬ 
rameter A = ^c(rr)/(rr)^. The single-particle Bogoli- 

ubov quasienergies E^(v n ) and Efc(v n ) are E^(v n ) = 
if — iv n — tk an d E^(y n ) = eg — iv n — t^. Next, fk and fk 
are band-renormalized hopping amplitudes fk = 2fgs7k 
and tk = 2fgs7k, where g b is the 3D bandwidth renor¬ 
malization factor 


g B 


,-i[ip(TT)-<p{r'r)\ 


|r— r'\—d 


(33) 


and 7k is the 3D lattice dispersion 7 k = cos {k x d x ) + 
cos(kyd y ) + cos {k z d z ), with d a (a = x,y,z), being the 
components of the lattice spacing vector d = r — r' with 
r and r', being nearest neighbors site positions. For the 
simple cubic lattice we have d a = a = 1. 

The quasiparticle energies If and eg are of Hartree- 
type and they are defined in the theory by relation 
e x = e x — g + Uriy + i (</(rr)), where y means orbital, 
opposite to x. The El low-temperature phase is charac¬ 
terized by the local excitonic order parameter A. With¬ 
out any loss of generality, we can suppose the case of the 
El state, with uniform real gap parameter A = A. The 
El state develops from local on-site electron-hole correla¬ 
tions. The expectation value, given in the expression of 
local El order parameter, could be calculated in the frame 
of path-integral method, as well as the fermion density 
averages of the respective band levels n x = (x(rr)x(rr)). 
We get a set of the coupled self-consistent equations for 
the El order parameter A, single-particle fermion densi¬ 
ties n x , and El chemical potential /i 


^ £[/(*£) +/TO] = 1 . 

k 


k 

A = V 

N ^ 

k 


/TO - /(E k ) 

V^+4A* 
/(£+) - f(E~) 
\/£k + 4A 2 


(34) 

(35) 

(36) 


Here /(e) = 1/ (e^ 6 + l) is the Fermi - Dirac distribution 
function, £k = —fk + eg + fk — £f is the quasiparticle 
dispersion and the energy parameters E£ and E ^ are 
defined as 


E 


± 

k 


1 

2 


tk + es — tk + i 



(37) 


In fact, the difference between Eqs. (PI) - ® and 
the Hartree-Fock results^ is in the presence of 


bandwidth renormalization factor gs- The calcula¬ 
tion of the factor gs (r — r') could be done effec¬ 
tively within the self-consistent-harmonic-approximation 
(SCHA) method-— AS In this approximation the quan¬ 
tum rotor description is reduced to classical Hamil¬ 
tonian one by the Feynmann-Kleinert minimization 
procedure - 47 ! 48 Our results for SCHA show that the fac¬ 
tor gB is equal identically to 1 at T = 0 as in the two- 
dimensional (2D) case (see in Ref. 133 for details). For 
higher temperatures, it differs from unity, but not much. 

The difference between the energy parameters in 
Eg. (1371) defines the charge transfer gap in the system 
A c = E£ — E^ = \/?k + 4A 2 (see in Ref. H3 for de¬ 
tails). 

The numerical solution of the system of equation 
Eas. (l3-H) - (15511 is discussed in details in Ref.([51), where 
finite-difference approximation method is used in nu¬ 
merical evaluations. Different values of t hopping am¬ 
plitude are considered there. The results are coincid¬ 
ing well with the previous Hartree - Fock (HF), im¬ 
proved slave-boson, and 2D constrained path Monte 
Carlo investigations j^ r 20 ■ 31 1 32 This good correspondence 
is ascribed to a rather weak band renormalization at 
T = 0 (in our case gs = 1). At finite temperatures, the 
particle number fluctuations are important and the band 
renormalization becomes necessary, especially, when ap¬ 
proaching from the band insulator (BI) high-temperature 
side. Indeed, as the numerical evaluations show, the tran¬ 
sition temperature Tei of the e-h pair formation is not 
vanishing for the case of the vanishing narrow band hop¬ 
ping t = 0. 

The exact solutions for the chemical potential could 
be obtained from Eas. (l34l) - (l36l) . both at the boundary 
of El transition (i.e., when A (Tei,U) = 0) and in the 
El state region ( A ( T < Tei,U) ^ 0). The results are 
discussed in Ref. (15 11). 

Meanwhile, it is also shown in Ref. Q that the ex¬ 
citonic BEC transition critical temperature T c is much 
smaller than the critical temperature Tei of the excitonic 
pair formation, in good agreement with previous theoreti¬ 
cal predictions - 15 ' 23-25 ! 47 ! 52 The self-consistent numerical 
solutions for Tei and T c are shown in Figs. [2] for t = —0.3. 
It is clear that, for the case of intermediate and strong 
interaction limits, the coherent BEC transition critical 
temperature T c is smaller of about two orders of magni¬ 
tude than temperature Tei ■ Contrary, in the very small 
interaction case, we have the coincidence of both tran¬ 
sition temperatures. This result is expected also from 
general considerations, because at small U/t we are in 
the BCS limit, which means that the pairing and con¬ 
densation occur simultaneously. 

A similar effect is found recently in Ref. 43, considering 
2D excitonic systems, where the exciton-superfluid tran¬ 
sition critical temperature is found much smaller than the 
critical temperature Tei of exciton pair formation. Such 
a reduction of the transition temperature, due to the co¬ 
herent pairing scheme, is given also in Refl52[ where the 
BCS-Bose crossover is studied in the 2D attractive Hub- 
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FIG. 2: (Color online) El transition critical temperature 
Tsi/t, (upper curve) and excitonic BEC transition critical 
temperature T c (lower curve). The data for the lower curve 
was taken from Refl5ll. 


bard model. Especially, a HF pair formation tempera¬ 
ture is estimated, which corresponds to the regime of the 
incoherent or local pairs (for a comparison see the cor¬ 
respondence to our Tei, and order parameter A, which 
is also local), in difference with superconducting pairing 
temperature, at which the coherent Cooper pair s start to 
be formed. The general idea used in Ref. {51) is based 
on the non local n.n. excitonic exchange correlation 
mechanism and Bogoliubov’s mean-field self-consistency 
assumption^ for the effective phase action obtained af¬ 
ter fermion integration procedure in the action in Eo. (127l) 
(this is represented by the right-lowest root in Fig. |l] 
given in the Section III Bl) . The obtained Hamiltonian for 
the phase sector is very similar to the classical Hamilto¬ 
nian one,^!~i£i£i with an effective phase-stiffness param¬ 
eter J that emerges after the fermion Wick averaging 
procedure^ The corresponding phase action is^i 



It appears that the non-zero value of quantity J is di¬ 
rectly related with the pairing gap A, since it is shown^i 
that J vanishes, when A = 0. Here, we present only the 
final analytical result for J (for details see in Ref. (51)) 


where z is the number of n.n. sites on the 3D lat¬ 
tice. The parameters Ai(k, k') and A 2 (k, k') entering in 
Eq.® are given by Ai(k,k') = (E+(k) -E^k'))” 1 ■ 

(E+(k)-E-(k')r\ A 2 (k,k') = (E-(k)-E-(k'))' 1 - 

(E _ (k) — E + (k')) 1 . Also, J is strictly positive for 
all the regions of the normalized Coulomb interaction 
parameter U/t. In addition, it follows from the an¬ 
alytical form of J, that the macroscopic phase coher¬ 
ence in the system is characterized by an energy scale 
J ~ t e ■ th/(t e + th) for all values of the Coulomb inter¬ 
action parameter U. This is related to the motion of the 
center of mass of e-li composed quasiparticle^ because 
t e ■ th/{t e + th) ~ (m e + TOft,) -1 . For the strong inter¬ 
action case, we are converging with the hard core Boson 
model, with the kinetic energy scale A t e ■ th/U (with A 
being the local excitonic order parameter). Thereby, it 
is shown in Ref. © that non-local correlations between 
the electrons and holes of different n.n. excitonic pairs, 
are related with the excitonic BEC condensation. Fur¬ 
thermore, in the frame of the quantum rotor model, the 
excitonic BEC transition probability function is derived 
and its temperature dependence is founds 

As a systematic continuation of theoretical study given 
in Ref. l(51j), we elaborate here on the analytical forms 
of normal, / and c -band, Green functions (incoherent 
and coherent) and coherent, anomalous excitonic Green 
function. Also, we will discuss in details, the spectral 
functions and density of states (DOS) corresponding, in 
the next. The numerical evaluations of calculated DOS 
functions are shown in the Section TlV Dl 


IV. SINGLE-PARTICLE DOS 

A. Normal and excitonic DOS functions 

The spectral density functions of the system of inter¬ 
acting exciton gas will determine the excitonic center- 
of-mass distribution related to the condensation in the 
low-temperature limit. Therefore the calculation of these 
functions represents an important task. Within our the¬ 
oretical approach we can access a variety of correlation 
functions in the system. We will concentrate now on the 
c and / -band normal spectral functions and also the ex¬ 
citonic anomalous spectral function, represented in terms 
of the initial operators. The c and / -band normal exci¬ 
tonic Green functions are 

G x , x (rr,rV) = )x(r / r')), 

(40) 


J = 


16A 2 ff ^ > 

z 2 N 2 ^ 


e(k)e(k') 

^+4A 2 


Ai(k, k') tanh 



—A 2 (k, k') tanh 



(39) 


the anomalous excitonic Green function is defined as 

G Cj f(rT,rV) = (c(rr)/(rV)). (41) 

After introducing the U(l) transformations, defined in 
Eq. (11141) , we will have the Green function’s decomposition 




















G x , x (rr,rV) = -(z(rr)l(rV)} • ( e -*Mrr)- v (rV)]^ 

(42) 

G c , f (rr,rV) = (g(rr)/(r'r')} • 

(43) 


The phase factors, appearing in the definitions of Green 
functions, define in fact a charge-bosonic propagator 
G z (yt,t't') as follows 

G 2 (rr,rV) = ( e -*M p G-v>( r V)]^ ( 44 ) 


Then, we pass to the Fourier space representation for the 
U(l) - transformed Green functions 


G XjX (rT,rV) = jjy Y ^x,x(k, v n ) x 

k.i/„ 

xe '[ k ( r -r')-‘'n(T-r')]^ 

G £ 4 (rr, rV) = j^Y ^£,f 1 "«) x 

k,I^n 

xe «[k(r-r')-j/ n (r-r')] 


(45) 


(46) 


and 


G 2 (rr,rV) 


1 

/3AT 


G z (k,w„)e i [ k(r - r ' ) -“’* (T - r ' ) ]. 

k,CcJ n 


(47) 


Furthermore, the Fourier transformation of the functions 
in Eas. (Hi?l) - (H51) will be rewritten as a convolution in the 
reciprocal k-space 


G X;X (k, Vn) 


1 

Jn 


Y^ G 2 (q, cc ra ) • G X)X (k 

q,oj n 


-q^n 


^n) 


(48) 


and 


G c ,f(k, v n )— qt\j 5 1 G z (q,o; n ) • G- f(k q, t'n w n ). 

q,u; n 

(49) 

It is worth to mention that frequency summations in 
Eas. (l48l) and (l49l) are over Bose-Matsubara frequencies 
uj n = 27rn//3. The Fermionic Green functions will be 

I 


calculated using the formalism discussed in the Sections 
nu and ED and also, functional derivation techniques.— 
Particularly, for the / and c -band Green functions 
G XiX (k, v n ), we get 


Gj, x (k, iv n ) = (i(k, v n )x(k, v n )) = 

= --. (50) 

. E-{v n )El{v n )- 

In general, the experimental observation of hybridization 
between the valence band and conduction band could be 
done by examining the ARPES spectra, which measure 
the spectral intensities just above and below the temper¬ 
ature Tei of excitonic pair formation. In ARPES ex¬ 
periments, one observes the imaginary part of the real 
retarded Green function, therefore the calculation of it 
represents a remarkable importance. Indeed, the single¬ 
particle density of states is related with the imaginary 
part of the retarded Green functions, thus we need to 
calculate real retarded function, which corresponds to 
the normal Matsubara Green function G^,x (k, iv n ). This 
could be done by the analytical continuation into the 
upper-half complex semi-plane ( v n > 0) of frequency 
modes iv n 


G^j(k, v) = G x ,s (k ,iv n ) \ iVn ^ v+il) . (51) 

The single-particle DOS is defined then as 

Pic.x (k, v) = - - ImG^ x (k, v) = (eg - f k - v) 2 x 

x<5 [(v 2 + A k v + .B k ) • (eg - ik - v)] , (52) 


where 


Ak — tk T Ik (-X eg (53) 


and 


Bk = ejeg + 4 tt [cos(fc x ) + cos (k y ) + cos(fc 2 )] 2 — 

—2 tlx [cos(fc x ) + cos (ky) + cos(fc 2 )] — 

—2teg [cos(fc x ) + cos (k y ) + cos(fc 2 )] — A 2 . (54) 

k-summed DOS will be 

Pic ,x {») = jj Y ( k ’ ^ ' ( 55 ) 

k 

The summations over the wave vectors in Eo. (15511 can 
be simplified by introducing the appropriate DOS func¬ 
tion for the 3D cubic lattice p 3 d{x) = jj 5Dk ~ 7k)- 
Then, it is not difficult to show that 


Pic,a (") 


r+3.0 


dxp 3D {x) 


[eg ~ t{ x ) - v] 2 
\/£ 2 {x) + 4 A 2 


S[v — E + (x )] 
\ey-t(x) -E+(x )| 


8[v — E (x)] \ 

\ey-t{x)-E~(x)\) 


(56) 


-3.0 
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Here, t(x) = 2 tx and t(x) = 2 tx, and the energy pa¬ 
rameters E±(x), are continuous versions of parameters 
defined in Ea. (l37ll . The density of state P 3 d{%), for the 
simple cubic 3D lattice, is given by 


The retarded function, which corresponds to it, is then^ 
Gfj( k, v) = G-j (k, iv n ) \ iVn ^ v+ir] . (59) 


P3d(x) 



(57) 


where 0(x) is the Heaviside step function, and K(x) 
is the elliptic function of the first kind: K(x) = 
JqG dtl / \/l — x 2 sin 2 t. 

For the anomalous excitonic Green function, we have 


G £ j(k, w n ) = (c(k, i/„)/(k, v n )) = 

A 

“ E*{y n )Ei{y n )- A 2 ' 


(58) 


J 


The single-particle DOS is forthcoming then as 


PcJ ( k > v ) 


= — — Im G'? -(k,;/) = 

7r 

A • S (v 2 + A\^u + . 


(60) 


Obviously, it has more simple form than the function in 
Ea. (l5i?l) . The k-summed DOS for excitons will be 


f P3D [Ai(t/)] p3D [A 2 (^)] | 

l IXi [ArM] | |X2 [A 2 (^)] | J ’ 


(61) 


where the dimensionless parameters Ai 2 {v) are given by 
following expressions 


A iM = 


A 2 (v) = 


(t + P) v — [e^t + e^tj + (t — t) v + [egt — eftj + 4tt|A| 


4 tt 


(t + t) v — + e ft j — (t — t) v + [ept — CftJ +4tt|A| 


4 tt 


(62) 


(63) 


and the functions Xi [Ai(i/)] (® = 1,2), in the denomina¬ 
tors in the right-hand side in Ea. flGTl) are 


Xi [Ai(i/)] = 2 (t + t) v + 8ttAi(v) - 2 (egt + ejtj 

(64) 


and 

X 2 [A 2 (i/)] = 2 (t +1) v + 8ttA 2 (v) - 2 (igt + cft) ■ 

(65) 

Now, turning to the convolution forms for total 
fermionic and excitonic Green functions in Eas. (148l) and 
(1491) , we need an explicit expression for the phase-bosonic 
Green function G z (k, ui n ). We will calculate it in the for¬ 
malism of the effective phase action given in the quantum 
rotor model, discussed earlier in Ref. [HTl, where we have 
derived the effective phase-only action £ e ff[y>] by inte¬ 
grating the fermions in the partition function in Eq. (|27D . 
In the following, we cast <5'eff , [¥ 5 ] into the quantum rotor 
representation^ 


B. The phase-stiffness DOS 

To proceed, we replace the phase degrees of freedom 
with the complex, unimodular field z(rr) = e ltp ^ TT \ 
which satisfies the time-periodic boundary condition 
z(rj3) = z(r0). The spherical constraint, imposed on a 
set of the unimodular variables is 1/N |z(rr)| 2 = 1. 

Now, we introduce these new variables z(rr) into the par¬ 
tition function in Ea. (127l) . in a way, consistent with the 
Faddeev-Popov ghost-field method^ 



xtf - e^ (rT) ) S[z- = 1 . (66) 

The phase-phase propagator G z (rr, y't') will be rewrit¬ 
ten in terms of z(rr)-variables as follows 

G z (rr,rV) = (z(rr)f(rV)). 

(67) 

Furthermore, the variables z(rr) play the role of the 
pliase-flux attached to the fermions (see discussions in 
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Section hd. In general case, the local expression of 
the phase-phase correlation function in Eg. (1771) is equal 
to unity, but, at very low temperatures (especially at 
T = 0), this law breaks down, because we have to con¬ 
sider the symmetry breaking related to the bosonic sec¬ 
tor, thus, critically, we have fluctuation form z(rr) = 
(e l¥, (rr)) + z(rr), and the unimodularity constraint for 
2 -field is violated. 

Indeed, in the very low temperature limit, considering 
the BEC of excitons, we have the spontaneous break¬ 
ing of local U(l) gauge-symmetry, related to the phase 
field, leading to the non-vanishing expectation value of 
the (e* v (rr)). In order to demonstrate this, we separate 
the single particle states k = 0 via the Bogoliubov dis¬ 
placement operation (see for details in Refsfll[5ll and (55) . 
This is so-called Bogoliubov phase coherence mechanism 
discussed in details in Ref. 1. Then, we write for the 
complex variables z(k,u> n ) 

z(k,uj n ) = PNi/joSk.oS^fi + z(k,w„)(l - <5 k ,o) x 

x (1o), (68) 


where Zq is the partition function of the non-interacting 
Bose sector 

+oo 

Z 0 = £ e-^( m "^) 2 . (73) 

m =—oo 

The summations, in Eos. (172l) and (1731) . run over topo¬ 
logical winding numbers m of the group U(l). The re¬ 
tarded bosonic Green function^ is related to the Mat- 
subara Green function, by the analytical continuation 

G?{k,u) = G z (k, iw n )|i w „->u+ir) ■ (74) 

And the k-summed DOS for bosons reads as 

P*{u) = --^ Im Gf(k,w). (75) 

k 

After non difficult algebraic manipulations and replacing 
the summation in Ea. (174l) by integration with the help 
of 3D density of states p-$D (x) = X)k $( x ~ 7 k) we get 


where i/jq is the BEC transition amplitude ipo = 
(z(k,w n )). Next, z(k,uj n ) are the excitation part£& (on- 
condensate) of effective bose-field. The Fourier trans¬ 
formation of the phase-phase propagator G z (rr, r'r') in 

Eg. (1771) is 

G z (kuj n ) = E {z(k,u n )z{k,uj n )). 

(69) 

We consider the expectation value (^(k, w n )f(k, uin)) in 
the local limit, i.e., when d = r — r' = 0 and t — t' = 0 
and we should draw the condensate part, by applying the 
transformation in Eg. (1681) . Hence, we have 

G z (k,uj n ) = — (z(k,u n )z{k,uj n )) 

= /3N\ip 0 \ 2 • <5 k ,o4„,o + G z ( k,w n ). (70) 


5[u- Kl{x)} 

^Jp 2 + AUJ (3 - x ) + 

S [u - k 2 (x)] 
v'/i 2 + 4D J (3 - x) ’ 

(76) 

where Ki{x) i = 1,2 are given by the following relations 
ki, 2 (x) = —fi ± \J\x 2 + 4 UJ (3 — x) and the stiffness pa¬ 
rameter J is given in Ea. (j38]l in Section Hill As it could be 
expected the Bosonic DOS function Eo. (176l) is negative 
p z (ui) < 0. This is consistent with the general consider¬ 
ations of the weakly non-ideal Bose gas— 

C. Spectral density functions and fermionic DOS 


AzM = - 


f+ °VA eE! 


Thereby, in Eq. (1701) we have defined the coherent macro¬ 
scopic state for the excitonic system in the low temper¬ 
ature limit, and an excitonic BEC is expected in the 
next. The Fourier Green function G z (k, ui n ), defined in 
Eg. (1771) . could be calculated within the quantum rotor 
model, and we give here only the final result (for details, 
see in Ref. l5lh . 


Furthermore, we separate the condensate modes 
(q = 0, ui n = 0} in Eos. (1751) and (1701) . We have 

Gx,x(k, Vn) — l^of * Gx ; x(k, V n ) T 

T njy 'y j G z (q, ui n ) ■ Gi j i(k — q, u n — ui n ) (77) 

q#0 

UJ n 7^0 


G z (k,ui n ) 


1 

7 _ 1 K) - 4Je(k) - A’ 


(71) 


where 7 ~ 1 (ui n ) is the inverse of the Fourier transforma¬ 
tion of the two-point phase-phase correlation function.— 
We have 


7 K) 


+ OO 

E 






( 72 ) 


and 

G c ,f(k, u n ) = |V>o| 2 • G 5 f(k, u n ) + 

+ ojy y j G z (<i,ui n ) • G~ f (k q, v n w n ). (78) 

q^O 

UJ n 7^0 


As we see, the normal and excitonic propagators are com¬ 
posed of two parts, one, responsible for the condensate 
state and the other- the on-condensate excitation part 
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(see discussion in RefjHH for the case of the pure Bose- 
gas). Note also, that first terms in the right-hand sides 
in Eqs. JZ3 and (1751) consist of the condensate-transition 
probability function I'i/’ol 2 ) multiplied with the fermionic 
propagators G X;X (k, v n ) and G ££ (k, v n ). 

Now, we are ready to calculate the analytical forms of 
the normal excitonic spectral functions A xx (k, u) (x = 
/, c) and anomalous excitonic spectral function A c j (k, v) 
and, later on, the profiles of the respective DOS, includ¬ 
ing states of the condensate. We introduce here the spec¬ 
tral functions A xx (k, v) and A c .f(k, v) that carries the 
same physical information as the correlation functions 
G XjX (k, v n ) and G Cjf (k, v n ). We have 


where n(e) = 1/ (e@ e — l) is the Bose-Einstein distribu¬ 
tion function. The proof of these relations is given in the 
Appendix ??. From the spectral functions we can obtain 
the corresponding DOS, by summing over the recipro¬ 
cal wave vectors k, hence, the total DOS are p x ,x(^) = 
7 ? Ek ^x,x(k, v) and p Cji (v) = ± E k A c,f ( k , v)- 
Furthermore, using the expressions for A XjX (k, v) and 
A C; f(k, v) in Eqs. (1571) and (1551) . we get for the total DOS 
functions 

PxA v ) = l^o| 2 • Px,x(^) + Px.xM (86) 

and 


„ , /■ +00 J .. / A XiX (k,^) 

— / dv . . (79) 

J-oo W n - v' 

G c f (k, v n ) = (80) 

The integration here, is over continuous frequencies. 
Note, that G XjX (k, v n ) and G c j (k, v n ) are total fermionic 
Green functions, including also the convolutions with 
bosonic parts. In the same way, we can introduce the 
spectral functions A z ( k, v), A XiX (k, v) and A~ £ (k, v ), as¬ 
sociated with the charge and the pure fermionic parts 
(without bosonic sector). They correspond respectively 
to the correlation functions G z (k,ui n ), G XiX (k, v n ) and 
G ££ (k, v n ). We have the following equations for these 
counterparts 

G a (k,o; w )= [ +0 ° dv , Az{KU '} , (81) 

J -OO *W n - V 1 

G x x (k, v n ) = (82) 

J-OO W n - v' 

and 


PcA”) = IV’ol 2 • P~ c j{y) + p z j(v), (87) 

where p^stiy) and p~ ^(y) are DOS, corresponding to the 
excitation part of the system and are given, as convolu¬ 
tions, in terms of continuous frequency modes (see Ap¬ 
pendix ??) 

/ -too 

dv'p z {y')Pi,z(y - v') 

-oo 

X [n (i/) + / {v - y')\ (88) 

and 

/ +oo 

dv'p z {v')p~ c i (v - v') 

-OO 

X [n (i/) + f{y- u')\ . (89) 

A key feature of the results in Eas. (l86l) and (1571) . is that 
we have separated DOS contributions coming from the 
condensate and excitation parts. We define also the total 
DOS function as 



du 


f Aq,{ ( k , 

iv n - v' 


(83) 


p( v ) = p^ v )- ( 9 °) 

x=f,c 


Using these definitions we get for the total spectral den¬ 
sity functions A XjX (k, u) and A c y(k, v) see Appendix ?? 


A x , x (k, v) = IV’ol 2 • A JjX (k, v) - 

~j7 5Z [ du ' A z ( q ’ ^)Ax,x(k - q, v - v') x 

x [n{y') + f(v - v')\ 

(84) 


and 


A c ,f(k, v) = |^ 0 | 2 ■ A z j{ k,i/) - 

[ dp'A z (q,u')A.^(k-- u') x 

x [n(y') + f(v - v ')\, 

(85) 


In the Section TlV D1 we present the numerical evaluations 
of all discussed DOS functions. 


D. Total DOS functions 

Employing Eqs. mm and (USD, we can obtain the 
explicite analytical expressions for DOS functions in 
Eas.(l86l) and (1571) . For the normal / and c -band ex¬ 
citonic DOS functions, we obtain 


Px.xM = |'0o| 2 • Px,x(l / ) 

P3d(x) 


r+3 


-u 


dx- 


J -s 4:\/p 2 + AUJ (3 — x) 

(Px.ic {v - Ki (x)) ■ [n {ki(x)) + f (v - Ki(a:))] + 

Px,x {y ~ «2 (x)) ■ [n (k 2 (x)) + / (v - k 2 (x))]} . (91) 
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For the anomalous excitonic DOS function, we have 

Pc,f{y) = \M 2 ■Pc,f( u ) 

_U [ + P^ x ) x 

7-3 V/i 2 + 4£/J (3 - a;) 

{pc,f - K 1 (®)) • [«(«iW) + / O' - «i(aO)] + 

Pc,t ( v ~ k 2 0*0) • [«(« 2 (*)) + f{v~ « 2 (a:))]} ■ (92) 

The numerical evaluations of calculated DOS functions at 
T = 0 are given in Figs. [3]- 0for t = —0.3. The presence 
of singularities in the integration region causes that we 
used an adaptive 21-point integration routine combined 
with the Wynn e-algorithm,— to calculate those integrals 
numerically. The accuracy for adaptive evaluations is 
achieved with an absolute error of order of ICR 4 and with 
a relative error of order of HR 7 . 

Particularly, in Figs. [3] and Fig. 0 we have presented 
purely fermionic normal single-particle (incoherent) DOS 
Px,x(w). We examine their behavior over the entire BCS- 
BEC crossover region (i.e., for different values of the 
Coulomb interaction U). An artificial Lorentzian broad¬ 
ening r] = 0.01 is used in numerical evaluations for the 
incoherent partial and total DOS functions of / and c - 
orbitals. The chemical potential values are inputed along 
the upper-bound Umax , where they are maximal (see in 
RefO 11) . The principal reason of it is that the BEC tran¬ 
sition amplitude tpo has no physical solutions along the 
lower-bound f.i m in of the chemical potential. On the other 
hand, the values n max are most convenient, because they 
are minimalizing the Hamiltonian of the system. Fur¬ 
thermore, as the reference for the coherent BEC tran¬ 
sition amplitude |'i/’o| 2 , we considered the self-consistent 
calculation-results from the work in Refl5p. where this 
function is calculated both analytically and numerically 
for different values of the Coulomb interaction parameter 
U and for different temperatures, including the zero tem¬ 
perature limit. We use here the numerical data, which 
are evaluated there. 

We see in Fig. 0 that at the small interaction limit, 
when 2 ^ U R 6 (i.e., the BCS limit), the incoherent 
excitonic DOS exhibits a BCS-like double-peak structure 
(see the panels I - III in Fig. [3]) and the peaks are sepa¬ 
rated with a well defined hybridization-gap. The princi¬ 
pal reason of it is the non-vanishing Hartree-gap A h ^ 0 
in the single-particle energy-spectrum discussed above, 
in the Section mu The hybridization-gap is proportional 
to the parameter U and it is increasing with the increase 
of parameter U. We observe also that the peaks becomes 
more separated when increasing of U. In the strong inter¬ 
action limit this displacement is stabilizing and we have 
practically constant value of the hybridization-gap, when 
further increasing the interaction (see the panels I - III, 
in Fig. 0. The results in Figs. [3] and 0are very simi¬ 
lar with the previous theoretical resultsj 15 ’ 20 ' 56 ’ 57 Espe¬ 
cially, they are close to those presented in RefsHH and 
|20| . where the partial incoherent / and c -band normal 


DOS functions and total DOS is calculated using HF and 
SO(2)-invariant slave boson approaches. 

In Fig. [5] (see the panels a - c), we have shown the co¬ 
herent condensate part of the anomalous excitonic DOS 
given by the first term in the right-hand side in Eq. (IA.4I) , 
corresponding to the fundamental state (k = 0). Dif¬ 
ferent values of the Coulomb interaction are considered, 
including the small and strong interaction cases. Here, 
we observe again the double-peak fermionic structure, 
but there is no the hybridization-gap for the small and 
medium values of the Coulomb interaction parameter 
(see the panel-a and panel-b in Fig. 0 , and we have fi¬ 
nite number of states for all values of the frequency modes 
v. This is due to the coherence effects and the presence 
of the coherent excitonic condensate at the fundamental 
mode v = 0. For higher values of U in Fig. [5] (see the 
panel-c), this double-peak structure in DOS is smooth¬ 
ing, but, in contrast to the incoherent normal DOS be¬ 
havior (see in Figs. [3] and 0]), here a small Mott-gap ap¬ 
pears at the very high values of the Coulomb interaction 
parameter (see the DOS curves for U = 8 and U = 9.6 
in the panel-c in Fig. 0. This is due to the fact, that 
the very strong Coulomb interaction has a destructive 
role on the condensate state and, in the large-R limit of 
interaction, we have the destruction of the excitonic con¬ 
densate, and a very small, Mott-type hybridization-gap 
is enhanced. In this region of the interaction, we have 
the coherent exciton DOS separation into two separate 
parts, similar to the case of the incoherent DOS functions 
in Figs. [3] and 0 

In the positive (negative) frequency regions, the DOS 
spectrum in Figs. 0 0 and 0 is slightly displacing 

and broadening into the direction of higher (smaller) fre¬ 
quencies, for both - normal (incoherent) single-particle 
fermionic, and coherent excitonic DOS functions and, for 
both, we observe also a gradual decreases in the DOS 
amplitudes, across the whole BCS-BEC crossover region, 
when increasing the Coulomb interaction parameter U. 
In Fig. 0 we have shown the total coherent fermionic 
DOS functions p(^), given by the Ea. fliJUl) . for different 
values of the Coulomb interaction parameter U and in 
the limit of zero temperature. An artificial Lorentzian 
broadening 77 = 0.01 is again used during the numerical 
calculations. Contrary the case of the incoherent DOS 
in Figs. 0 and 0 the coherent total fermionic DOS in 
Fig. 0 shows a different, gapless behavior. Here, again, 
as in the case of the incoherent DOS functions, for small 
and intermediate interactions U, we have typical double¬ 
peak fermionic structure in the DOS, which is smoothing 
and disappearing totally in the strong interaction lim- 
ite. But, in difference with incoherent DOS functions, 
here we have not the presence of the hybridization-gap 
in the spectra and we have always have a finite number 
of states for all values of the frequency modes (see in 
Fig. 0. The reason of this gapless DOS behavior could 
be related to the strong coherence effects between two 
bands, which is due to the presence of the phase stiffness 
mechanism considered here. It is worth to mention that 
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FIG. 3: (Color online) Single-particle normal DOS functions 
(incoherent) for different values of the Coulomb interaction 
parameter U and for the case T = 0. In panels I-III, x = f 
and x = c DOS structures are shown, and total DOS is plotted 
for t = —0.3. 


FIG. 4: (Color online) Single-particle normal DOS functions 
(incoherent) for different values of the Coulomb interaction 
parameter U and for the case T = 0. In panels I-III, x = f 
and x = c DOS structures are shown, and total DOS is plotted 
for t = —0.3. 
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FIG. 5: (Color online) Condensate DOS function |?/>o| 2 -p g f f° r 
different values of the Coulomb interaction parameter U and 
for the case T = 0. The case t = —0.3 is considered here, and 
the values of the functions \tpo\ 2 were taken from the work in 
Ref. U 



V 


FIG. 6: (Color online) Total phase coherent DOS function 
p(v) given in Eg. (1901) for different values of the Coulomb inter¬ 
action parameter U (U = 4, U = 6, U — 8, U = 9, U == 10.4 
in the figure) and for the case T = 0. The case t = —0.3 
is considered here, and the values of the functions |V'o| 2 were 
taken from the work in Ref. l5ll . 


another gapless-type behavior in the DOS spectrum is 
found recently in Refj60j, where this effect is associated 
with metallic charge-density-wave phase and is driven by 
strong electron correlations. In Fig. [Gj all figures are 
combined together, in the way, to see the total DOS evo¬ 
lution with variation of the interaction parameter U. We 
see clearly, how the tuning of the interaction parame¬ 
ter U , affects the general DOS behavior in the model, 
by reducing the DOS amplitudes with increasing of U 
and reducing also the number of states at the Fermi level 
p(v = e F ) when increasing U. Thereby, in the case of 
normal fermionic DOS (both incoherent and coherent), 
presented in Figs. [3] [I] [5] and Fig. 0 the double-peak 
structure disappears for U > 6, signaling the appearance 
of the SM (BEC) limit of the transition. 

In Fig. 0 a comparison is given for the incoherent and 
coherent c -band normal DOS behaviors. Three differ¬ 
ent values of the Coulomb interaction U are chosen (see 
the panels -a, b and c). We see clearly in Fig. 0 how 
the coherence effects are reducing the incoherent DOS 
amplitudes and DOS spectra become broader, along fre¬ 
quency axis, and also, there is no gap near the Fermi level 
{p(v ~ e F ) ± 0). 

In Figs. [8] and 0 we have presented the tempera¬ 
ture dependence of the c -band normal fermionic (Fig. 0) 
and anomalous excitonic DOS functions (Figs.0). They 
are given by the first terms in the right-hand side in 
Eas. dA.4l) and (l92l) . In the left-panel in Fig. 0 we 
have presented the temperature dependence of the single¬ 
particle DOS I'i/’ol 2 Pc,c(v), which corresponds to the fun¬ 
damental state k = 0, and for U = 6. The case t = —0.3 
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is considered. The corresponding values of the amplitude 
ipQ of BEC transition, are taken again from the work in 
RefIFlI. In the right-panel in Fig. [9l the same function 
is plotted for the case U = 9. As we see in Fig. [51 the 
hybridization-gap is still open for all values of the tem¬ 
perature. So the system is always an insulator. 

In Fig. El the temperature dependence of the phase- 
coherent anomalous excitonic condensate part of the 
DOS is presented (i.e., k = 0, and without excitation 
part) for U = 2 and U = 6. We observe in all Figs. [5] 
and [H that the temperature has a destructive effect on 
the DOS amplitudes. 

We realize also, from Figs. [5] and E that the anoma¬ 
lous excitonic DOS functions vanish at the temperatures, 
that are far away from the region of the El transition of 
about two orders of magnitude (compare the tempera¬ 
ture scales in Figs. [Hand [9] with those given in Fig. [2] 
in the Section ED). This result, could be regarded as a 
good proof of the theory elaborated in ReflU, and we can 
theoretically clearly state that the excitonic BEC and El 
states are not the same phases of matter. We see also, in 
Figs. [H an d El that the normal single-particle DOS per¬ 
sists for a rather large values of temperature, than the 
anomalous condensate DOS, because of the presence of 
the hybridization-gap in the low-energy spectra. 


V. FINAL REMARKS AND CONCLUSIONS 

We have studied 3D system of conduction band elec¬ 
trons and valence band holes in the frame of the ex¬ 
tended Falicov-Kimball model. We have implemented 
the path-integral formalism, in which the Coulomb in¬ 
teraction term is expressed in terms of U(l) quantum 
phase variables ip conjugated to the local particle num¬ 
ber, providing a useful interpretation of the problem. In 
Section ED we have shown that at low temperatures, the 
electron-hole system becomes unstable with respect to 
the formation of the excitons at T = Tei , and the local 
gap A is present in the excitation spectrum, controlled by 
the Coulomb interaction parameter U/t, which gives the 
relevant energy scales for the excitonic insulator state. 
Here, as a result of the spontaneous symmetry breaking, 
an expectation value of (e l¥> ) ^ 0 appears, which is signal¬ 
ing of the presence of the phase coherence in the system. 
Furthermore, pairing and condensation are not generally 
the same, as it was admitted in the literature, except the 
weak interaction limit, when we have a BCS-likc con¬ 
densation of excitonic pairs. However, in the excitonic 
system with the strong pairing, we have the situation, 
where the pairs are strongly bound, but are uncorrelated 
one with each other, until they become phase coherent 
at temperatures T < T c . This situation was studied in 
details in ReflEI. 

We have evaluated the normal and anomalous exci¬ 
tonic spectral functions for the / and c -band. We have 
determined, both analytically and numerically, density 
of states (DOS) spectra, governed by the pure-fermionic 
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FIG. 7: (Color online) Incoherent and coherent normal DOS 
functions for the c -band and for different values of the pa¬ 
rameter U (U = 4 in a, U = 6 in b and U = 8 in c). The case 
t = —0.3 is considered here. 
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FIG. 8: (Color online) Temperature dependence of the normal 
c -band DOS function for the case t = —0.3 and for different 
values of the Coulomb interaction parameter U. 




FIG. 9: (Color online) Temperature dependence of the 
anomalous excitonic DOS function for the case t = —0.3 and 
for different values of the Coulomb interaction parameter U. 


part (due to the condensate modes k = 0), and the ex¬ 
citation spectra. We have shown that there is a usual 
hybridization-gap in the normal (incoherent) / and c - 
band DOS structures. Contrary, in the case of the coher¬ 
ent normal DOS functions (presented in Fig. [6]) this gap 
is lacking, and there is always a finite number of states 
at all frequency modes v. We associate this result to the 
strong coherence effects, that are present in the system 
at low temperatures. 

In the anomalous excitonic DOS structure we have 
found that the hybridization-gap is absent for the weak 
and intermediate values of the Coulomb interaction pa¬ 
rameter U and this is due to the presence of the coherent 
excitonic condensate and strong coherence effects. A very 
small gap is opening in the spectra in the strong interac¬ 
tion limit, signaling the destruction of the coherence and 
condensate state. 

The excitonic phase coherence may be evidenced by the 
coherence of their light emission, which can be studied 
by interferometry measurements ^ Therefore, measure¬ 
ments of the intensity of the line-shape of the excitons 
decay (by emitting photons, upon electron-hole recombi¬ 


nation) may be a powerful probe of the DOS spectra in 
the excitonic systems^ 

However, a final remark that should be featured is 
such, that for the experimental determination of the ex¬ 
citonic Bose-Einstein condensate, the ARPES measure¬ 
ments should be provided at temperatures much lower 
than temperatures at which the excitonic insulator state 
is being deternrinedi^This is important for the achieve¬ 
ment of macroscopic phase coherece between excitonic 
pairs, and the strong coherence effects are manifesting at 
the very low temperatures. 

As a continuation of our studies, we would like to con¬ 
sider the role of the charge-density-wave-like excitations 
in the excitonic systems, within the EFKM model and 
consider the temperature effects. It is especially interest¬ 
ing to find out how the coherent excitonic DOS will be 
affected in the case of presence of such elementary exci¬ 
tations. This will give us a more complete picture on the 
excitonic phase transition scenario. 
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Appendix: A convolution for DOS 


Here, we give a short derivation of Eqs. (l84l) and (l85l) 
presented in the Section II V Cl To this end, we use the 
denitions in Eas. lRiTl) - (1551) and the convolution form of 
the Green functions given in Eqs. m and CS) . Then, 
it is easy to see that the following identity holds for the 
integrals of the normal a; = /, c -band fermionic spectral 
functions. 


d y Ax,x (k, v) =| ^ o|2 


dv A *’* (k,Z/) + 


+ 


PN 


E 

q^0,a;ri.^0 

l 


dv / dv' 


iui n — iv n + u' — V 


A z (q, v') • Ajj (k — q, v - v'). 

(A.l) 


Furthermore, for calculating the Matsubara sums over 
bosonic frequencies co n , we will use the property of the 
Bose-Einstein distribution function n (e) (see the Section 
IIV Cl) of complex argument 

n (iv n + v — v') = —f {v — v '), (A.2) 

where / (e) is the Fermi-Dirac distribution function. In 
Ea. dA.2l) . we used the fact that v n are even Fermion- 
Matsubara frequencies u n = ^(2n + l). Next, we sum 
the Bosonic Matsubara frequencies in Ea. dA.il) amd we 
rewrite the equality in Ea. dA.il) in the following form 
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A x ,x (k, v) 


l^oI 2 • 4r,x (k, v) 


[ dv ' Az ^ q ’ v> ) ^’ 5 ( k - q » v - v ') • \ n ( v ') + f(y- v ')\ ■ 

q#0 J 


(A.3) 


The derivation of the convolution form, given in Eq. (18511 , 
is exactly the same. 

Summing over the wave vectors k in Ea. (IA.3|) . we will 
have the forms of the normal / and c -bands density 
of states p X}X (v ) and also the excitonic DOS function 
p c j (is) given in Eas. (15Sl) and (1571) in the Section HV Cl 


1. Normalization conditions for DOS 


We have a composed nature of the interacting electron, 
given in the Section UD thus, we have various sum rules, 
corresponding to different counterparts of the total coher¬ 
ent DOS functions. As a general rule, the partial normal 
DOS, and anomalous excitonic DOS functions p X:X (is) 
and p c j (is) satisfy following normalization conditions 


r*+oo 


dvpif (is) = U{, 


/• + OO 


disp c ^ c (is) = n c , 


(A.4) 


and 


/ +oo 

dispef (is) = 1 

-OO 


(A.5) 


For the total normal band DOS function p(is), given in 
Eg. (1501) . we have 


f* + 00 


disp (is) = 1 


(A.6) 


The functions pay(is) and p~ j( is ) in Eas. (1551) and (1551) are 
excitation parts^ of the respective total DOS functions 
in Eos. (1561) and (1571) . i.e., without condensate modes k = 
0, therefore, the following normalization conditions hold 


/ -t-oo 

dvp~ t ~ t (v) = rif - \tp 0 \ 2 , 

-OO 

/ +oo 

dispzfi(is) =n c - IV’ol 2 (A.7) 

-OO 


and 



1 - IV'ol 2 ■ 


(A.8) 


In contrary, the integral over the charge-DOS function 
p z (u), given in Eg. (1751) . vanishes, because it is antisym¬ 
metric: p z (u>) = —p z (—ui). Then 


/• + 00 


dujp z (ui) = 0 


(A.9) 


The norms of the full spectral density functions in 
Eos. (l84l) and (l85l) depend on the constraints in the charge 
and spin bosonic sectors. 
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